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ABSTRACT 

In A Bayesian Approach to Locating the Red Giant Branch Tip Magnitude (PART I),' a new technique 
was introduced for obtaining distances using the TRGB standard candle. Here we describe a useful com- 
plement to the technique with the potential to further reduce the uncertainty in our distance measurements 
by incorporating a matched-filter weighting scheme into the model likelihood calculations. In this scheme, 
stars are weighted according to their probability of being true object members. We then re-test our modified 
algorithm using random-realization artificial data to verify the validity of the generated posterior probability 
distributions (PPDs) and proceed to apply the algorithm to the satellite system of M31, culminating in a 3D 
view of the system. Further to the distributions thus obtained, we apply a satellite-specific prior on the satel- 
lite distances to weight the resulting distance posterior distributions, based on the halo density profile. Thus 
in a single publication, using a single method, a comprehensive coverage of the distances to the companion 
galaxies of M3 1 is presented, encompassing the dwarf spheroidals Andromedas I - III, V, IX-XXVII and XXX 
along with NGC147, NGC 185, M33 and M31 itself. Of these, the distances to Andromeda XXIV - XXVII 
and Andromeda XXX have never before been derived using the TRGB. Object distances are determined from 
high-resolution tip magnitude posterior distributions generated using the Markov Chain Monte Carlo (MCMC) 
technique and associated sampling of these distributions to take into account uncertainties in foreground ex- 
tinction and the absolute magnitude of the TRGB as well as photometric errors. The distance PPDs obtained 
for each object both with, and without the aforementioned prior are made available to the reader in tabular 
form. The large object coverage takes advantage of the unprecedented size and photometric depth of the Pan- 
Andromeda Archaeological Survey (PAndAS). Finally, a preliminary investigation into the satellite density 
distribution within the halo is made using the obtained distance distributions. For simplicity, this investigation 
assumes a single power law for the density as a function of radius, with the slope of this power law examined 
for several subsets of the entire satellite sample. 

Subject headings: galaxies: general — Local Group — galaxies: stellar content 



1. INTRODUCTION 

The Tip of the Red Giant Branch (TRGB) is a well estab- 
lished standard candle for ascertaining distances to extended, 
metal poor structures containing a sufficient red giant popu- 
lation. Its near constant luminosit y across applicable stellar 
mass and metallicity ranges (see llben & Renz ini 1983) arises 
due to the prevailing core conditions of these medium-mass 
stars as core helium fusion ensues. Their cores lack the nec- 
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essary pressure to ignite immediate helium fusion on the de- 
pletion of their hydrogen fuel and so they continue to fuse 
hydrogen in a shell around an inert, helium ash core. This 
core is supported by electron degeneracy, and grows in mass 
as more helium ash is deposited by the surrounding layer of 
hydrogen fusion. On reaching a critical mass, core helium 
fusion ignites, and the star undergoes the helium flash before 
fading from its position at the tip of the Red Giant Branch, to 
begin life as a Horizontal Branch star. Due to the very similar 
core properties of the stars at this point, their energy output is 
almost independent of their total mass, resulting in a distinct 
edge to the RGB in the Color-Magnitude Diagram (CMD) of 
any significant red giant population. 

With the TRGB standard candle applicable wherever there 
is an RGB population, it is an obvious choice for obtain- 
ing distances to the more sparsely populated objects in the 
Local Group and other nearby groups where Cepheid Vari- 
ables seldom reside. Even when Cepheids are available, the 
TRGB often remains a more desirable alternative, requir- 
ing only one epoch of observation, and facilitating multiple 
distance measurements across an extended structure. Good 
agreement between TRGB obtained distances, and those ob- 
tained using Cepheid Variables as well as the much fainter 
RRLyrae Variables have been confirmed by Sala ris & Cassisil 
d 19971). with discrepancies of no mo re than ~ 5% (see also 
Tammann, Sandage, & Reindll 120101 for an extensive list of 
distance comparisons utilizing the three standard candles). Of 



2 



Conn et al. 



the satellites of M31, many are very faint and poorly popu- 
lated and thus have poorly constrained distances which propa- 
gate on into related measurements concerning the structure of 
the halo system. Hence, a technique for refining the distances 
that can be applied universally to all halo objects, whilst ac- 
curately conveying the associated distance errors, has been a 
long sought goal. 

In 'A Bayesian Approac h to Locating the R ed Giant Branch 
Tip Magnitude (PART I)' - IConn et al .1(1201 ll) , hereafter Paper 
I, we reviewed the challenges of identifying the TRGB given 
the contamination to the pure RGB luminosity function (LF) 
typically encountered. We also outlined some of the meth- 
ods that have been devised to meet these challenges since the 
earlies t approach, put forward by iLee. Freedman. & M adore 
(1993). We then introduced our own unique Bayesian Ap- 
proach, incorporating MCMC fitting of the LFs. This ap- 
proach was essentially the base algorithm, designed to easily 
incorporate priors to suit the task at hand. Here we present 
the results of an adaptation of that algorithm, intended for use 
on small, compact objects - specifically the dwarf spheroidal 
companions of M3 1 . Once again, we ut ilize the data of the 
Pan- A ndromeda Archaeological Survey dMcConnachie et aTl 
2009), a two-color (/' = 770 nm, g' — 487 nm) panoramic 
survey of the entire region around M31 and M33 undertaken 
using the Canada-France-Hawaii Telescope (CFHT). The tip 
is measured in i' band where dependance on metallicity is 
minimal. Following a recap of the base method in §|2] we 
int rodu ce the a fore mentioned new adaptations to the method 
in {33. H and in $3.2l we describe the results of tests intended to 
characterize the modified algorithms performance as well as 
check the accuracy of its outputs. In addition, ^3.31 outlines 
t he a pplication of a further prior on the satellite distances. 
§14. II presents the results of applyi ng th e modified algorithm 
to the companions of M31, while ^4. 21 details the m etho d by 
which the object-to-M3 1 distances are obtained and jj4.3l uses 
the obtained distances to analyze the density profile of these 
objects within the halo. Conclusions follow in §0 

2. A RECAP OF THE BASE METHOD 

In Paper I, we introduced our 'base' method, whereby the 
LF of a target field was modeled by a single, truncated power 
law (the RGB of the object of interest) added to a represen- 
tative background polynomial. The location of the truncation 
(the TRGB) and the slope of the power law were set as free 
parameters of the model, with the best fit derived using an 
MCMC algorithm. The functional form of the background 
component was modeled by directly fitting a polynomial to 
the LF of an appropriate background field, and then scal- 
ing the polynomial to reflect the expected number of back- 
ground stars in the target field. The resulting model was 
then convolved with a Gaussian of width increasing in pro- 
portion to the photometric error as a function of magnitude. 
The posterior distribution in the tip magnitude returned by the 
MCMC, which thus already incorporates the photometric er- 
ror, is then sampled together with Gaussian distributions rep- 
resenting the distribution in the absolute magnitude of the tip 
(MJRcb _ _3 44 + o.05) and the distribution in the extinc- 
tion (A A ± 0.1 A A ) to give a final posterior distribution in the 
distance. The mode of this distribution is then adopted as the 
distance to the object, with the ±l<x error calculated from the 
portion of the distribution lying on the far and near side of the 
mode respectively. 

A more detailed discussion of the assumptions and ratio- 
nale behind the base method is provided in paper I , but the 



reader should again be made aware of the most fundamen- 
tal assumptions it entails. At the heart of the calculations of 
course is the choice of the absolute magnitude of the tip and 
its associated uncertainty. We adopt the values of this param- 
eter stated abo ve based on the value derived for the SDSS 
i-band in iBellazzinil (120081) . noting the near-identical band- 
pa ss characterist ics of the MegaCam i-band filter as detailed 
by iGwvnl (I2010I) . We adopt som e what smaller uncertainties 
than those derived by Bellazzini (2008) following the same 
argument as McConnach ie et all ((2004) that the quoted uncer- 
tainty in the absolute magnitude of the tip is conservative, and 
it is a systematic error effecting all distance measurements in 
an identical way. As almost all applications of the distances to 
the satellites are concerned with their relative positions to one 
another and M31, this component of the error is of minimal 
importance. Nevertheless, it often forms the major component 
of the quoted errors in our distances. 

Mention should also be made as to the effects of metallic- 
ity and internal reddening within the objects under study as 
well as the zero-point uncertainty in the PAndAS photometry. 
Whilst there is a metallicity dependence of Mj RGB , (though 
minimal when compared with other bands), it is only really 
an issue for mor e metal-rich targets (e.g. [Fe/H] > -1, see 
IBellazzinil 2008 Fig. 6) and thus will primarily effect mea- 
surements to the large, diverse systems such as M3 1 itself and 
M33. But the TRGB for more metal rich populations is fainter 
than that for their metal-poor counterparts and thus it is this 
metal-poor population component which dominates the mea- 
surement. A similar situation is encountered with the internal 
reddening present in the objects under study, where the vast 
majority of objects, chiefly the dwarf spheroidal galaxies, are 
almost completely devoid of such effects. Those objects most 
strongly effected are the large, well populated systems which 
will provide ample signal from the least effected stars on the 
near side of the system, for a good distance determination. 
The uncertainty in the zero-point of the photometry is consis- 
tent throughout t he survey at approximately 0.02 magnitudes 
(llbata etalJl2012l) . 

Lastly, a brief discussion of the distance posterior distribu- 
tions themselves is warranted. As noted above, they are pro- 
duced by the sampling of the distribution of possible tip posi- 
tions (as generated by the MCMC and with photometric errors 
incorporated) along with sampling of the Gaussian distribu- 
tions representing the uncertainties in the foreground extinc- 
tion (A A ) and in the absolute magnitude of the tip {Mj RGB ). 
Specifically, 500, 000 possible distances are drawn to form the 
distance PPD, where for each draw k, the distance modulus /j 
is: 



A a {k)- 



mJ rgb (k). 



(1) 



where each of mJ RGB (K), A a {k) and MJ rgb (k) are the val- 
ues drawn from the uncertainty distributions in the tip posi- 
tion, foreground extinction and absolute magnitude of the tip 
respectively. The foreground extinction and its uncertainty 
varies from object to object but the error in the absolute mag- 
nitude of the tip is a systematic error as already discussed. In 
using this method, there are two situations that can be encoun- 
tered. The first is that the object is very well populated and 
the tip position is thus well constrained with a narrow PPD. 
In such instances, the uncertainty in Mj RGB far outweighs any 
other contributions to the error budget and is almost solely re- 
sponsible for the width of the distance PPD. In the second 
situation, the object is poorly populated and the tip position 
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PPD is very wide and typically asymmetric. If the LF pop- 
ulation is not extremely low, the uncertainty in Mj RGB will 
contribute noticeably to the distance PPD, otherwise the dis- 
tance PPD will essentially depend solely on the uncertainty 
in the determined tip positions. Hence whilst some of the 
smaller contributions to the distance uncertainties are omitted 
from the calculations, their over all effects will be washed out 
by the contributions from these two principle sources of error. 

3. ADDITION OF A MATCHED FILTER 

3.1. Matched Filtering using Radial Density Profiles 

With the introduction of our method in Paper I, it was 
stressed that one of its greatest attributes was its adaptabil- 
ity to the prior knowledge available for the object of inter- 
est. When applying the method to compact satellites, there is 
one very conspicuous attribute that can be incorporated into 
the prior information constraining the model fit - namely, the 
object's density as a function of radius. The simplest way to 
achieve this is with the addition to the algorithm of a matched- 
filter weighting scheme, wherein the weighting is matched to 
the specific data by accounting for the data within the filter 
itself. 

The successes of iRockosi et al.l ((2002) using a matched fil- 
ter in colour-magnitude space to identify member stars of 
globular cluster Palomar 5 amidst the stellar background pro- 
vide the inspiration for our technique. They make use of the 
characteristic RGB of the globular cluster to weight stars as 
to their likelihood of being cluster members. To achieve such 
a goal, a matched filter can be created by binning the CMD of 
the field in which the cluster lies into a 2D matrix and then di- 
viding that matrix by a similarly created background matrix. 
Stars found in the densest regions of the resulting matched 
filter CMD are then assigned the highest weight, being the 
most likely cluster members. In this way, they can greatly 
improve the signal-to-noise ratio (SNR) with respect to that 
of their original, unmodified data and are able to trace tidal 
streams from the globular cluster well into the surrounding 
background. Hence we have applied a similar approach to 
weight field stars fed to the MCMC in terms of their prob- 
ability of being object members. In our case however, the 
stars proximity to the object's center provides the basis for the 
weighting scheme, with the inner most stars being the most 
likely to be actual object members as opposed to background 
stars, and so a one-dimensional matched filter is sufficient. 

The first step in implementing our weighting scheme is to 
ascertain a model of stellar density as a function of radius 
specific to the object of in terest. For thi s purpo se, we employ 
the best fits presented in iMartin et al.l d2012h for the dwarf 
spheroidal satellites, wherein the optimal ellipticity e, position 
angle (PA), half-light radii (r;,) and object centers are given 
for exponential density profiles fitted to each satellite. For 
the two dwarf ellipticals, in the case of NGC147 we a ssume 
e = 0.44 and PA = 28° as specified bv lGeha etail (120101) and 
we derive the r„ manually as 10', which produces the best fit 
profile to the data when coupled with the other 2 parameters. 
For NGC185, we adopt e = 0-26 and PA = 41° based on 
the findings of Hodge (1963) and once again derive the r/, 
manually, this time as 6'. For both NGC147 and 185, we 
employ the object centers derived from the 2 Micron All Sky 
Survey (2MASS - ISkrutskie etail l2006). With the ellipticity, 
position angle, half-light radius and object center know, we 
can proceed to produce a weighting scheme proportional to 
the density profile p of the object, where p is of the form: 



p(r e ) = e"r (2) 

where R = is the scale radius and r e is the elliptical radius 
at which the star lies, as now defined. With the PA and object 
center of the object known, a rotation of coordinates is used 
to define each star's position (x',y') with respect to the center 
of the ellipse. The projected elliptical radius r e of the ellipse 
on which the star lies is then: 




where the y' axis is assumed as the major axis of the ellipse. 

Whilst Eq. [2]gives us the functional form of our weighting 
scheme, it is further necessary to define the absolute values 
of the weights given to each star, so as to scale them appro- 
priately with respect to the background density p/, g . This is 
achieved by insuring that the area under the function p(r € ) be- 
tween any imposed inner and outer radius limits is set equal 
to the number of signal stars in the observed region. Hence, 
our weighting scheme is ultimately defined by: 

W(r e ) = Se^r (4) 

with, 

_ (p total ~ Pbg)><A 

2nR(\ - e)[(e^)(R + r imer ) - (e^XR + r outer )] 

(5) 

where p tota i is the density of stars in the observed region be- 
fore subtraction of the background density, A is the area of the 
observed region which is either an ellipse in the (usual) case 
that Turner = or an elliptical annulus otherwise. r mner and 
r outer are the inner and outer cutoffs respectively of the range 
of r e values observed. 

In Fig. [1] the result of our fitting procedure as applied to 
the sparsely populated dwarf spheroidal Andromeda X is pre- 
sented. In this case, stars out to r e = 0.15° are fitted, with 
no inner cutoff radius imposed. Whilst most of the satellites 
are too poorly populated for blending to be an issue, in the 
case of several, the stellar density counts at the inner most 
radii drop off in spite of the predicted counts from the fitted 
density profile. This is a good indicator of blending or over- 
crowding in those radii which can hinder the accuracy of the 
photometry for the affected stars and so in such cases, these 
inner radii are omitted. This was the case with Andromeda 
HI (rtnner = 0.0175°), Andromeda V (r imer = 0.011°) and 
Andromeda XVI (r,„„<, r = 0.005°). For the dwarf ellipticals 
NGC147 and NGC185, it was found beneficial to avoid the 
inner regions all together, with the presence of a wider range 
of metallicities in these regions degrading the contrast of the 
RGB tip. Similarly, an outer cutoff radius was chosen for 
these objects inside of 3 r/, to help sharpen the tip disconti- 
nuity, so that for NGC147, r ilmer = 0.28° and r outer = 0.33° 
and forNGC185, r inner = 0.18° and r outer = 0.26°. M31 and 
M33 are treated similarly to the dwarf ellipticals but with still 
thinner annuli so that any w eigh ting is unnecessary. They are 
discussed in more detail in 84.11 

With regard to the actual likelihood calculations used at 
each iteration of the MCMC, these are undertaken not by sim- 
ple multiplication of the likelihood for each star by the re- 
spective weight, but by physically adjusting the relative pro- 
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Elliptical Radius (degrees) 

Fig. 1 . — Radial Density Profile (proportional to Object Membership Prob- 
ability) for Andromeda X. The error bars represent the Poisson Error in the 
density for each bin, with each bin representing an elliptical annulus at the 
stated radius. Hence the inner most annuli have the smallest areas and thus 
the largest error bars. Note that this binned density distribution is for com- 
parison only and has no bearing on the fit. The background level is marked 

'bg: 



portions of the RGB and background components of the lu- 
minosity function. Up until now, we have assumed a generic 
LF and calculated the likelihood contributions from each star 
from this single LF. But in reality, the outer regions of the field 
are more accurately represented by a shallow-signal/ high- 
background LF while the innermost stars obey a LF which 
has almost no background component. Hence using the radial 
density profile obtained above, we can essentially build an in- 
dividual LF for each star, tailored to suit its position within 
the object. In practice, this is achieved with almost no extra 
computational effort, as the background and signal can be nor- 
malized separately and only the signal component is changed 
by the MCMC at each iteration so that the background com- 
ponent need only be generated once. The two components are 
normalized to contain an area of unity and then the bin of each 
corresponding to the star's magnitude is scaled according to 
the ratios of the star's weight and the background level when 
its contribution to the model likelihood is calculated by the 
MCMC. 

The result of the incorporation of this extra prior informa- 
tion is a marked improvement in the performance of the algo- 
rithm for the more sparsely populated targets. In such objects, 
the RGB component is typically overwhelmed by non-system 
stars, even with the most carefully chosen field size. This can 
greatly diminish the prospects of obtaining a well constrained 
tip measurement. This is apparent from figures[2]and[3]which 
show the luminosity function and corresponding posterior dis- 
tributions before and after the application of the matched fil- 
ter to the dwarf spheroidal Andromeda X. With the matched 
filtering applied, the great majority of non-system stars are 
severely suppressed, revealing clearly the RGB component, 
which in turn provides much stronger constraints on the loca- 
tion of the tip, as evidenced by Fig. [3] Herein lies an example 
of the power of the Bayesian approach, where a single prior 
can cast the available data in a completely different light. 

3.2. A Test for the refined algorithm 




19.5 20 




Fig. 2. — Best fit model to the luminosity function of Andromeda X, ob- 
tained with the addition of matched filtering. The top figure shows the best 
fit overlaid on the un-modified LF (i.e. histogram created without the weight- 
ing afforded by the matched filter). The bottom figure shows the same best 
fit model after applying the weighting. A field radius of 0.15° was used to 
generate the LF histograms, wherein each star contributes between and 1 
'counts', depending on its proximity to the field centre and the density profile 
of the object. 

In §2.3 of Paper I, the results of a series of tests were pre- 
sented that characterized the performance of our original al- 
gorithm given a range of possible background density levels 
and LF populations. Here we present the results of similar 
tests applied to our new, matched-filter equipped algorithm, 
but with some important differences. Most fundamentally, the 
way our artificial test data is generated is quite different. As 
we are now concerned with the position of each star in the 
field, a distance from field centre must be generated for each 
star. To do this, we have randomly assigned a radial distance 
to each star, but weighted by a circularly symmetric (e = 0) 
exponential density profile . Further to this, the magnitudes of 
our stars are now generated directly from our convolved LF, 
so that photometric error as a function of stellar magnitude is 
incorporated. 

The other important change from the previous tests con- 
cerns the way in which the artificial luminosity functions 
are populated. Whereas in the former tests all of the sam- 
pled stars were drawn from the model LF within the one 
magnitude range 20 < m star < 21, in the current tests the 
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Fig. 3. — Posterior distributions obtained for Andromeda X before (top) and 
after (bottom) the application of the matched filter. For the 'before' case, 
a circular field of radius 0.05° (2.143 X r;,) has been chosen, specifically to 
provide the most possible signal with the least possible background contam- 
ination. For the 'after' case, the same LF as presented in Fig. [2]is used. 

stars are drawn from within the much larger magnitude range 
actually utilized for our satellite measurements, namely 
19.5 < m slar < 23.5. Hence a 100 star LF in these tests 
for example corresponds to a much smaller sample of stars 
than in the tests described in §2.3 of Paper I. Aside from 
these critical differences, the current tests are undertaken 
and presented as per the previous publication, with mea- 
surements of the average sigma and tip offset given for 
each combination of background level (/) .vs. number 
of stars (ndata) where / = 0.1, 0.2, 0.9 and ndata = 
10, 20, 50, 100, 200, 500, 1000, 2000, 5000, 10000, 20000. 
The results are presented in figures[4]and[5]respectiveiy. 

Examination of the figures reveals the expected trend of in- 
creased one-sigma error and tip offset with increasing back- 
ground height and decreasing LF population levels. Once 
again, there is very good agreement between the derived er- 
rors and the actual offsets obtained. Most importantly, it is 
clear by comparing these results with those of Paper I, that 
the Matched Filtering has greatly diminished the effects of 
the background contamination, as exemplified by the much 
gentler increase in \<j errors and offsets with increasing back- 
ground star proportion. 
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Fig. 4. — A grey-scale map of the one-sigma error in tip magnitude obtained 
for different combinations of background height and number of sources. The 
actual value recorded for the error (in kpc) is overlaid on each pixel in red. 
Each value is the average of twenty 50000 iteration runs for the given back- 
ground height/ LF population combination. 
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Fig. 5. — A grey-scale map of the offset of the measured tip value from the 
true tip value obtained for different combinations of background height and 
number of sources. The actual value recorded (in kpc) is overlaid on each 
pixel in red. Each value is the average of twenty 50000 iteration runs for the 
given background height/ LF population combination. 



3.3. An additional Prior 

In addition to our density matched filter, a further prior may 
be devised so as to constrain our distance Posterior Probability 
Distributions (PPDs) in accordance with our knowledge of the 
M3 1 halo dwarf density profile. The expected fall off in den- 
sity of subhaloes within an M3 1 sized galaxy halo is not well 
constrained. The largest particle simulation of an M31 sized 
dark m atter halo to date, the Aquarius Project (Springel et al. 
2008), favored the density of subhaloes to fall off following 
an Einasto Profile with r_2 = 200 kpc and a = 0.678, and 
furthermore identified no significant dependence of the rela- 
tionship on subhalo mass. For the specific cas e of the satellites 
within the M3 1 halo, Richards on et alj d201 lb found a relation 
of p oc r~ a where a - \ a. better fit to the data, drawing largely 
from the PAndAS survey, although this does not take into ac- 
count the slightly irregular distribution of the survey area. We 
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adopt this more gentle density fall off with radius giving us a 
more subtle prior on the satellite density distribution and note 
that a may be changed significantly without great effect on 
our measured distances. 

So in effect, we assume a spherical halo centered on M31, 
such that p(sat) oc r and integrate along a path through the 
halo at an angle corresponding to the angular displacement on 
the sky of the satellite from M3 1 . This yields an equation of 
the form: 

d 2 

P(d) <* . (6) 
V(t/ 2 + 779 2 - 2d x 779 x cos(0)) a 

where a = 1, 779 kpc is the distance to M31 and P(d) is 
the relative probability of the satellite lying at distance d (in 
kpc) given an angular separation of 6 degrees from M3 1 . Note 
that this produces a peak where the line of sight most closely 
approaches M31, and that P(d » 779) is approximately pro- 
portional to d. The equation is normalized between limits ap- 
propriate to the size of the halo. 

We thus generate a separate prior for the probability as a 
function of distance for each satellite, tailored to its specific 
position with respect to M31. The effect of the prior is to 
suppress unlikely peaks in the multi-peaked posterior distri- 
butions obtained for certain satellites, while leaving the peak 
positions unaffected. As such, the prior has very little effect 
on single peaked distributions, whatever the angular position 
and distance of the satellite it represents. The distance prior 
applied to the Andromeda XIII distance PPD is shown in Fig- 
ure |6]for illustration. 
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Fig. 6. — The Distance Prior applied to Andromeda XIII (solid line; a = 1). 
The distribution gives the likelihood of the satellite existing at a particular 
distance, given an angular separation on the sky of 8.5° from M31 (the halo 
centre), and assuming a distance of 779 kpc for M3 1 . The distribution peaks 
where the line of sight traverses the inner most region of the halo, and flattens 
out at large distances due to the increasing volume of the halo subtended by 
the unit of solid angle observed. The same prior with a = 2 is shown as a 
dashed line for ref erenc e. Whilst this value for alpha is in closer agreement 
with the results of j|4.3l we deliberately adopt the less restrictive a = 1 prior, 
so as not to suppress the probability of satellites in the outer halo too greatly. 



4. A NEW PERSPECTIVE ON THE COMPANIONS OF 

M31 

4. 1 . Galaxy Distances 



The PAndAS Survey provides us with a unique opportu- 
nity to apply a single method to a homogeneous data sample 
encompassing the entire M31 halo out to 150 kpc. The data 
encompasses many dwarf spheroidals, along with the dwarf 
ellipticals NGC147 and NGC185, and of course the M31 disk 
itself with additional fields bridging the gap out to the com- 
panion spiral galaxy M33, some 15° distant. Of these objects, 
the vast majority have metallicities [Fe/H] < -1, so that any 
variati on in the absolut e magnitude of the tip is slight. In- 
deed, iBellazzinTI ([2008) suggests that for such metallicities, 
the variation in the region of the spectrum admitted by the 
CFHT i' filter is perhaps less than in Cousins' I. Perhaps of 
greatest concern are the cases of M31 and M33, which will 
contain substructure at a variety of metallicities. In this case 
however, the more metal rich portions will exhibit a fainter 
TRGB than those in the regime [Fe/H] < -1, such that the 
brightest RGB stars will fall within this regime. 

In this section we present distance measurements to these 
many halo objects, culminating in Fig. [10l a three- 
dimensional map of the satellite distribution, and Table [2j 
which presents the satellite data pertinent to our distance mea- 
surements. Figures [TTI and [121 present the distance posterior 
distributions obtained for every object in this study. It has 
been common practice in the majority of TRGB measure- 
ments to quote simply the most likely distance and estimated 
lcr uncertainties, but this throws away much of the informa- 
tion, except in the rare case that the distance distribution is 
actually a perfect Gaussian. On account of this, as well as 
providing the actual distance PPDs themselves for visual ref- 
erence, we also provide the same information in condensed 
tabular form, where the object distance is given at 1% incre- 
ments of the PPD, both for the prior-inclusive cases (as in 
Figures QT| and [T2l > and for the case in which no prior is in- 
voked on the halo density. Note that for M3 1 , no halo density 
prior is applied and so this column is set to zero. A sample of 
this information, as provided for Andromeda I, is presented 
in Table [1] The reader may then sample from these distribu- 
tions directly rather than use the single quoted best fit value, 
thus taking into account the true uncertainties in the measure- 
ments. 

Due to the large number of objects studied, it is not practical 
to discuss each in detail within this publication. For this rea- 
son, Andromeda I will be discussed in further detail below as 
a representative example, followed by two of the more prob- 
lematic cases for completeness. Firstly however, we describe 
the exceptional cases of M31 itself and M33. 

M31 and M33, due to their large extent on the sky and the 
variety of substructure in their disks require a slightly differ- 
ent approach to that used for the other objects in this study. As 
was the case for NGC147 and NGC185, it was necessary to 
define a thin elliptical annulus so as to limit as much as possi- 
ble the amount of substructure from other radii contaminating 
the LF. For both M31 and M33 such a thin annulus was used 
that any weighting with respect to the elliptical radius of the 
stars was trivial and so no weighting was used. For M31, an 
ellipticity of 0.68 was adopted, with PA = 37°. The inner and 
outer elliptical cutoff radii were set to 2.45° and 2.5° respec- 
tively. To check for any inconsistencies in the TRGB location 
across the whole annulus, it was divided up into NE, NW, SE 
and SW quarters and then the distance measured from each 
quarter, giving distances of 782+}^, 782+J|, 775^° and 781^ 
kpc respectively. It is tempting to associate the slightly lower 
distance to the SE quadrant with the effects on the LF of the 
Giant Stellar Stream, though the distance is still within close 
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TABLE 1 

Tabulated Distance Posterior Distribution: 



Percentage 


Distcin.ce (kpc - no density prior) 


Distsnce (kpc) 


i 


684 


687 


2 


688 


692 


3 


691 


695 


4 


693 


697 


5 


695 


699 


6 


697 


701 


7 


698 


702 


8 


699 


703 


9 


700 


704 


10 


701 


705 


100 


820 


820 



Distance posterior probability distributions for Andromeda I given at 1 % 
intervals for the case of no halo dens ity prior (column 2) and with the 
angle-specific prior outlined in j|3.3l applied (column 3). 



agreement with the other 3 quadrants, such that all 4 are per- 
fectly consistent. Hence, the distance was remeasured using 
the whole annulus to give 779 + } | kpc. This is in good agree- 
ment both with the findings of McConnac hie et ail (|2005A) 
(785^) utilizing the TRGB, and the more recent determi- 
nation bv lRiess. Fliri. & V alls-Gabaud (2012) using Cepheid 
Variables (765^). 

For M33, we employ an ellipticity of 0.4 as used by 
iMcConnachie et aTl (I2005AI) . but find a position angle of 
PA = 17° in closest agreement with the data. Inner and 
outer elliptical radii of r mner = 0.75° and r outer - 0.9° were 
adopted to give a very sharp discontinuity at the location of 
the tip. After applying an appropriate color-cut, the qualifying 
stars were fed into our algorithm to give a distance of 820^ 

kpc. This distan ce is in good agreement with that of 809^ 
kpc obtained by McConnachie et al. (2005 A) and yields an 
M33 to M31 distance of 214^ kpc. It is interesting to note 
that a variety of quite different M33 distances exist in the 
literature, with d e rived distance moduli ranging from 24.32 
dBrunfhaler et al.l (120051) - 73 kpc - water masers) through 
24.92 dBonanos et all (120061) - 964 kpc - Detached Eclips- 
ing Binaries). Indeed, the variety of standard candles utilized 
would suggest that M33 provides an ideal e nvironment for 
calibrat ing the relative offsets between them. IMcConnachie! 
(2005B) suggests that the dispersion of M33 distances in the 
literature is tied to an inadequate understanding of the extinc- 
tion in the region of M33. Most measurements, including 
those p resented here, use the Galactic extinc tion values de- 
rived by Schlegel, Finkbeiner, & Davis ( 1998), although these 
do not account for extinction within M33 itself and are cal- 
culated via an interpolation of the extinction values for the 
surrounding region. Nevertheless, the elliptical annulus em- 
ployed in our approach will act to smooth out the field-to-field 
variation that might exist between smaller regional fields. 

4.1.1. Andromeda I - Example of an ideal luminosity function 



It would seem prudent to illustrate the performance of our 
new method by presenting the results for a range of the dwarf 
spheroidals from the most populated to the least populated. 
Hence Andromeda I, the first discovered and one of the two 
most highly populated of these objects is the obvious place 
to start. The field employed for our Andromeda I distance 
measurement incorporated stars at elliptical radii between 
0° < r £ < 0.3° and, after removal of stars outside of the range 
19.5 < io < 23.5 and beyond our chosen color-cut, yielded 
a star count of 4375. The CMD for this field is presented in 
Fig. |7]a. This figure color-codes the stars in the CMD as per 
the color distribution in the inset field and plots them so that 
those innermost within the field (and hence those accorded 
the highest weight) are represented by the largest dots. In the 
case of Andromeda I, the RGB is so dominant over the back- 
ground that our density matched filter is hardly necessary and 
hence does little to improve the already stark contrast. It is not 
surprising therefore that the distance and uncertainty obtained 
are almost identical to those obtained by the base method as 
presented in Paper I. Andromeda I is thus confirmed at a dis- 
tance of 727+ [7, which allows us to derive a similarly accurate 
separation distance from M31 of 68^^ kpc. 

4.1.2. Andromeda XV - Example of a multi-peaked distance PPD 

As an example of a dwarf spheroidal of intermediate size, 
we present the comparatively compact Andromeda XV. Far 
from being the tidiest example of the many intermediate sized 
objects covered in this study, Andromeda XV provides some- 
thing of a challenge. Examination of Fig. [8]reveals a gradual 
rise in star counts when scanning from the top of the CMD 
color-cut faintward toward the Andromeda XV RGB and a 
correspondingly broad range in the possible tip locations in 
the tip magnitude PPD. Indeed, two peaks are prominent in 
the distance PPD of Fig. [8] c, with the distribution mode at 
626 kpc (our adopted distance) and the lcr credibility inter- 
val spanning fro m 591 kpc to 705 kpc as a consequence of 
the second peak. Ilbata et all ([2007) determine this object to 
lie at a distance of 630+^ kpc, which would correspond to 
a tip magnitude of approximately mJ RCB = 20.56 assuming 
Mf RGB = -3.44. This is in excellent agreement with the 
m JRGB _ 20.57+J}^ recovered by this study. iLetarte et all 
(2009) however derive a distance of 770*22 kpc which places 
it toward the far edge of our 99% credibility interval on the 
distance (see Fig. [8] c). This measurement was derived 
after 3 stars that had been found to lie close to the An- 
dromeda XV RGB tip in the former investigation were identi- 
fied as Galactic foreground stars, following measurements of 
their radial velocities obtained with the Deep Imaging Multi- 
Object Spectrograph on Keck II. Of these stars however, none 
lie within 2' from our object center, by which point the max- 
imum possible weighting has already dropped to below 10%, 
meaning that even the highest weighted of these 3 stars will 
have minimal effect on the likelihood calculation. This would 
then suggest that each of these 3 stars have magnitudes con- 
sistent with belonging to the Andromeda XV RGB. 

4.1.3. Andromeda XIII - Example of a very poorly populated 
luminosity function 

Andromeda XIII is among the most sparsely populated ob- 
jects targeted by the current study and it is important to re- 
alize that it is impossible to obtain distances to such objects 
with small uncertainties using the TRGB standard candle, un- 
less of course one of the few member stars can be positively 
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identified as being right on the brink of core Helium fusion. 
Nevertheless, though large uncertainties are inevitable, an ac- 
curate estimation of those uncertainties is still achievable, and 
this is the aspiration of the method here presented. Distances 
to Andromeda XI and XIII have been obtained with higher 
accuracy using RR Lyrae stars as a Standar d Candle with pho- 
tomet ry from the Hubble Space Telescope ( Yang & Saraiedini 
2012). In the case of Andromeda XI, the tip magnitude iden- 
tified by our method agrees well with the distance identified 
by that study, but in the case of Andromeda XIII, a brighter 
star in the central regions of the field causes some confusion. 
Indeed in such a sparsely populated field it is quite difficult to 
apply any effective density-based weighting scheme. Never- 
theless, after sampling the tip magnitude PPD (Fig. |9](b)), to- 
gether with those for the absolute magnitude of the tip and the 
extinction in this region of sky to obtain a sampled distance 
PPD, and multiplying that distribution with the angle-specific 
halo density prior as is standard for all our measurements, we 
are able to produce a dist ance PPD (Fig. |9](c ) ) in g ood agree- 
ment with the findings of lYang & S araiedinil (120121) . 

4.2. Determining the distances from M31 

Once a satellite's distance from Earth is determined, it is 
straightforward to determine the distance from M3 1 using the 
cosine rule: 

r = ^d 1 + (c/m3i) 2 - 2ddM3icos(8) (7) 

where r is the satellite's distance from M31, d is the distance 
of the satellite from Earth, d*f3i is the distance of M31 from 
Earth and 9 is the angle on the sky between M3 1 and the satel- 
lite. For convenience, we use a small angle approximation 
equating 6 with its M31 tangent plane projection and note 
that any displacement of r is insignificant due to the size of 
the lcr errors. If the uncertainty in distance to both M31 and 
the satellite takes on a Gaussian distribution, it is straightfor- 
ward to determine the error in the satellite-M3 1 separation by 
adding the individual errors in quadrature. While it is rea- 
sonable to approximate the M31 distance uncertainty distri- 
bution as a Gaussian, the same cannot be said for each of the 
companion satellites. Hence once again it is more appropri- 
ate to sample values from the individual distance probability 
distributions. Thus, a histogram of r values for the satellite is 
built up by sampling d and o?m3i from their respective distri- 
butions over many iterations. This brings to the fore an im- 
portant consideration: there is an integrable singularity in the 
resulting distribution at the closest approach distance to M3 1 
(r c = dMiisin(d)) as shown below. 

The probability distribution for the Satellite-to-Earth dis- 
tance P(d) is related to that of the Satellite-to-M31 distance 
P(r) as follows: 

Sd 

P(r) = — P(d) (8) 
or 

From Eq. [7] and further noting that the Satellite-to-Earth dis- 
tance corresponding to r c is d c = dM3\cos(8) we have: 



5r d - d c 
which allows us to derive: 

P(r) = r P(d) (10) 



thus producing the singularity at r — r c . In practice, after 
factoring in the Gaussian distribution in dmi, this results in 
a sharp peak at the minimum possible Satellite-to-M31 dis- 
tance when dealing with the more asymmetric Satellite-to- 
Earth distance probability distributions. Hence when consid- 
ering the distribution of satellites as a function of distance 
from M31, one can either take the distances as determined 
directly from Eq. fusing solely the most likely distance from 
the Satellite-to-Earth distance distributions, or the whole dis- 
tance probability distribution for a satellite can be allowed to 
influence the calculations, as accomplished via sampling. The 
final result can be quite different, depending on the choice. 

4.3. A First Approximation of the Satellite Density Profile 
within the Halo 

In the completed PAndAS Survey, we have for the first time 
a comprehensive coverage of a galaxy halo, with a uniform 
photometric depth sufficient to identify even the compara- 
tively faint satellite companions. In addition, in this paper 
we have provided distances to every one of these objects, all 
obtained via the same method. We are thus presented with 
an excellent opportunity to study the density of satellites as a 
function of radius within a Milky-Way-like halo. 

As hinted at in the previous section, obtaining an accurate 
picture of the Satellite Density Profile (SDP) is not a trivial 
task. The first major consideration is to devise a way of fac- 
toring in the selection function. Comprehensive though the 
survey coverage is, it is not symmetric, and it is not infi- 
nite. Secondly, the choice of model for the SDP is not ar- 
bitrary. Whether a simple, unbroken power law is sufficient 
is not immediately clear. Furthermore, does it even make any 
sense to treat the halo as a radially symmetric, isotropic dis- 
tribution? A glance at the obvious asymmetry in Fig. [T0l(a) 
would suggest otherwise. Nevertheless, for a first approxima- 
tion it is reasonable to consider what the best fitting radially- 
symmetric, unbroken power law to the SDP would be. 

The PAndAS Survey covers approximately 400 square de- 
grees of sky and is roughly symmetric about the center of the 
M31 disk but with a major protrusion in the South East to 
encompass the M33 environs. For the purpose of obtaining 
an accurate measure of the survey coverage of the halo as 
a function of radius, as well as factoring in the actual sur- 
vey borders, an inner ellipse was also subtracted where the 
presence of the M31 disk has made satellite detection more 
difficult. Both the outer survey borders and the inner cut-off 
ellipse are plotted in Fig. |T0j The inner cut-off ellipse has an 
eccentricity e = V0.84 and is inclined with the semi-major 
axis angled 51.9° with respect to the x-axis (77 = 0). The 
dwarf galaxies M32 and Ml 10 lie inside this ellipse as do the 
somewhat dubious satell ite identifications Andromeda VII I 
and Andromeda IV (see lFerguson. G allagher. & Wyse 2000), 
hence their omission from the data presented in Table[2] With 
the inner and outer boundaries suitably delineated, the proce- 
dure then was to determine what fraction of halo volume at 
a given radius f(r) would fall within these boundaries once 
projected onto the M31 tangent plane. This was achieved by 
implementing the Even-Odd Rule on the projections of uni- 
formly populated halo shells. 

Having determined /(r), we can proceed to determine the 
required normalization for a power law of any given a, allow- 
ing us to use the power law directly as a probability distri- 
bution. Setting the problem out in terms of probabilities, we 
require to determine the probability of each tested M31-to- 
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20-8 20.82 20.84 20.86 20.88 20.9 20.92 19.5 20 20.5 21 21.5 22 22.5 23 23.5 



Proposed i Q tip magnitude L 

Fig. 7. — Andromeda I: a): Color-coded CMD representing the weight given to each star in the field. Only stars within the red selection box with magnitudes 
19.5 < i'o < 23.5 were fitted and hence color-coded. The second, fainter RGB lying toward the redder end of the CMD is that of the Giant Stellar Stream which 
passes behind our Andromeda I field. The inset at top right shows the field with the same color-coding and acts as a key. The field is divided into 20 radii bins 
following a linear decrease in density from the core (blue) to the field edge (purple). Stars marked as a purple 'x' lie outside of the outer elliptical cutoff radius 
r ouur . Stars marked as a black '+' are artificial stars used in the estimation of the background density and are ignored by the MCMC; b): Posterior Probability 
Distribution for the TRGB magnitude. The distribution is colour coded, with red indicating tip magnitudes within 68.2 % (Gaussian 1-sigma) on either side of 
the distribution mode, green those within 90 % and blue those within 99 %; c): Weighted luminosity function of satellite with superimposed best-fit model in red. 
A star at the very centre of the satellite contributes 1 count to the luminosity function while those further out are assigned some fraction of 1 count in proportion 
with the satellite's density profile 
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Fig. 8.- — Andromeda XV: a): Same as Fig. |JJ(a) but for Andromeda XV.; b): Same as Fig. |7](b) but for Andromeda XV; c): Sampled Distance Posterior 
Probability Distribution, obtained by calculating the distance 3 million times, each time randomly drawing on the tip magnitude, absolute magnitude of the tip 
and extinction from their respective probability distributions. The distribution is colour coded, with red indicating possible distances within 68.2 % (Gaussian 
1-sigma) on either side of the distribution mode, green those within 90 % and blue those within 99 %. Note that the large uncertainty in the absolute magnitude 
of the RGB tip is primarily responsible for the much smoother appearance of the distance PPD (c) compared with the tip PPD (b). 




Fig. 9. — Andromeda XIII: All figures as per Fig. [8] but for Andromeda XIII. The distance derived by Yang & Sarajedini (2012) is plotted in (c) along with 
error bars for comparison. 
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TABLE 2 
M31 Satellite Parameters: 

Distance and associated parameters of M31 and its companions. All distance measurements utilize the data from the Pan- Andromeda Archaeological Survey 
I McConnachie et al. 2009), and have been obtained using the method presented in this paper. A value of Mj RGB = -3. 44 ± 0.05 is assumed for the absolute 
magnitude of the RGB tip in CFHT MegaCam i-band, based on the v alue identified for the SDSS i-band ( Bellazzini 2008) and justified for use here by the color 
equations applicable to the new MegaCam i-band filter i Gwyn 2010). Values for the extinction in MegaCam ('-band have been adopted as = 2.086 X E(B - V) 
for the same reasons, with uncertainties taken as ±10%. Note that the uncertainties in the M31 distance are based on the sampled distributions while the quoted 
value is that derived directly from the earth-distance as per Eq. [7] The last column gives alternative distances from the literature. TRGB derived distances are 
quoted wherever possible. 



Source 


Distance Modulus 


E(B-V) 


Distance (kpc) 


M31 Distance (kpc) 


Literature Distance Values (kpc) 


M31 


n * *,-+no5 
24.46!°;°= 


0.062 


779+.;* 




785+— TRGB; McConnachie et al. (2005A) 
784+i 7 , RC; Stanek & Garnavich (1998) 
765+ 28 Ceph; Riess, Fliri, & Valls-Gabaud (2012) 


And I 


24 3 j +0.05 

ZH - -0.05 


0.054 


797+I8 
'-17 


68- 


731 + j? TRGB; Conn et al. (2011) 
735+?? TRGB; McConnachie et al. (2004) 


And II 


24 00+ 005 


0.062 


630! J 3 


195! 2 ° 


634+j' TRGB; Conn et al. (2011) 
645+ J' TRGB; McConnachie et al. (2004) 


And III 


24 W°- U3 
zqou -0.07 


0.057 


723+.^ 


86!- 


749+ 24 TRGB; McConnachie et al. (2005A) 


AndV 


74 3 =+0.06 
-0.07 


0.125 


742+H 


113!' 


774! 28 . TRGB; McConnachie et al. (2005A) 


And IX 




0.076 


600!^ 


182-66 


765! 24 TRGB; McConnachie et al. (2005A) 


AndX 


24.13+°;?* 


0.126 


670+2* 


130!» 


667 - 738 TRGB; Zucker et al. (2007) 


And XI 


24 41 +0.08 
z ^ 1 -0.32 


0.080 


763+ 2tJ 

""-106 


102!j 49 


740 - 955 TRGB; Martin et aL (2006) 
735+!?, RRLv; Yang & Saraiedini (2012) 


And XII 


24 84+ a(w 
z ^- -0.34 


0.111 


928+ 40 

?zo -136 


181- 


825 + _f 5g TRGB; (MCMC without MF) 
740 - 955 TRGB; Martin et al. (2006) 


And XIII 


24 40 

-0.49 


0.082 


760 

' -154 


1 15 2 


890+ j* TRGB; (MCMC without MF) 
740 - 955 TRGB; Martin et al. (2006) 

839+jS RRLv; Yang & Saraiedini (2012) 


And XIV 


n A c/i+0 06 

24.50_ 056 


0.060 


793 -179 


161™ 


630 - 850 1RGB; Maiewski et al. (2007) 


And XV 


23.98 !2 


0.046 


626 -35 


1 HA +46 

1 74 -32 


*;">r\+60 TTir "n. ti <- n „+ „1 i"1AA-7\ 

60 i *^* J ^' Ibata et al. (200/) 
770+™ TRGB; Letarte et al. (2009) 


And XVI 


23.39+°-!'; 

-0.14 


0.066 


476!^ 


319+^ 

-27 


525^° TRGB; Ibata et al. (2007) 
525!^° TRGB; Letarte et al. (2009) 


And XVII 


24 31+0-11 
zt.ji _ 0Qg 


0.075 




67 20 

"'-24 


794+ 40 . TRGB; Irwin et al. (2008) 


And XVIII 


25 42+ lIU7 


0.104 


1214!^ 


457 +39 


1355+S 8 . TRGB; McConnachie et al. (2008) 


And XIX 


24 = 7 +o.08 
"^•-"-0.43 


0.062 


821^8 


ns*l b 


933+IJ TRGB; McConnachie et al. (2008) 


And XX 


24 3 5+°- 12 


0.058 


741+42 


12%+f 


802! 2 ^ 7 TRGB; McConnachie et al. (2008) 


And XXI 


74 <;q+0.06 
'-0.07 


0.093 


827 -25 


135+ 8 


859+jj TRGB; Martin et al. (2009) 


And XXII (Tri I) 


24.82+°;° 7 


0.075 


920+ 32 

y -139 


275+ 8 


7g4 +239 T/RGB; Martin et al. (2009) 


And XXIII 


24.37+°;- 


0.066 


748— 


127! 7 


733+ 23 TRGB; Conn et al. (2011) 
767+ 44 HB; Richardson et al. (2011) 


And XXIV 


24 77+0-U7 

z +-"-o.io 


0.083 


898! 2S , 


169+ 2y 


600+S HB; Richardson et al. (2011) 


And XXV 


24 33+0.07 
"■"-0.21 


0.101 


736+- 

-69 


'"-10 


812+ 4 5 HB; Richardson et al. (201 1) 


And XXVI 


24 39+°- 55 


0.110 


754+218 
IJ -164 


103! 234 


762+ 42 HB; Richardson et al. (2011) 


And XXVII 


25 40+0.07 
■"•^-1.03 


0.080 


1255+« 4 


482!° 25 


827! 47 HB ; Richardson et al. (2011) 


And XXX 1 ' (Cass II) 


24 17+0-10 
Z4 - -0.26 


0.166 


681$ 


145!f 


565!- TRGB g-band; Irwin et al. (2012) 


NGC147 


24 26+°" U6 
z ^- zo -0.06 


0.173 


712- 


118!- 


675+ 27 , TRGB; McConnachie et al. (2005 A) 


NGC185 


23 96+ lIU7 
-0.06 


0.182 


620!- 


181- 


616+ 2 * TRGB; McConnachie et al. (2005A) 


M33 


24 57+°- 05 

ZH - J '-0.05 


0.042 


820! 2 ° 


210!* 


809+ 24 TRGB; McConnachie et al. (2005 A) 
964+^ DEB; Bonanos et al. (2006) 



Note : Extinction values are for object centers. Actual calculations apply individual corrections to each member star according to their coordinates. 

Distance derivation methods: TRGB = Tip of the Red Giant Branch; Ceph = Cepheid Period-Luminosity Relation; RR Ly = RR Lyrae Period-Luminosity 
Relation; RC = Red Clump; HB = Horizontal Branch; DEB = Detached Eclipsing Binary 

' Andromeda XXX is a new disc overy, and will als o be known as Cassiopeia II, being the second dwarf spheroidal satellite of M3 1 to be discovered in the 
constellation of Cassiopeia - see Irwin et al. ( 12012ft 
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Fig. 10. — Three views of the M31 Neighborhood: a) A view of the satellites of M31 along the y-z plane. The conic section illustrates the extent of volume 
covered by the PAndAS footprint as a function of distance from Earth; b) A view of the satellites of M3 1 in the x-y plane, revealing their true positions on the 
x-y plane after removing the effects of perspective (assuming the distances quoted in column 4 of Table [2). Note that Andromeda XXVII lies directly behind 
NGC147 in this plot and is not labeled.; c) A 3D view of the satellites of M31. The satellite positions on the PAndAS footprint are indicated (i.e. with perspective 
conserved) along with the z- vector giving distance from the M3 1 - centered tangent plane. The central ellipse indicates the approximate area of the survey where 
satellite detection is hindered by the M31 disk; Note: The perpendicular bars on relevant axes indicate 100 kpc intervals. 
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Fig. 11. — Distance Posterior Distributions for dwarf spheroidal satellites And I - III, And V and And IX - XIX. The distributions are color-coded with red, green 
and blue denoting 1-cr (68.2%), 90 % and 99 % credibility intervals respectively. The credibility intervals are measured from either side of the highest peak. 




Fig. 12. — Distance Posterior Distributions for dwarf spheroidal satellites And XX - XXVII and And XXX, dwarf elliptical satellites NGC147 andNGC185, and 
major galaxies M31 and M33. The distributions are color-coded with red, green and blue denoting l-cr (68.2%), 90 % and 99 % credibility intervals respectively. 
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object distance (henceforth simply 'radius') r given a power 
law with slope a: 



P(r\a) = 



(11) 



where k is the normalization constant and r,„„, < r < r max . 
Using the assumed spherical symmetry, we then have: 



f(r) 



P(r\a) I r 2 sin6ddd(f>dr = 1 (12) 

Jr„,i„ -JO JO 



so that 



4nf(r) ^ "'" k?- a dr = 1 
Hence, for a given radius at a given a, we have: 



k(r, a) 



Anf{r) 



Ji-a 



3 — a 



(13) 



(14) 



The calculation of the likelihood for a power law of a given 
slope a may be simplified by noting that for any given radius, 
f(r) and hence k act to scale the probability in an identical 
way whatever the value of a. Thus the dependance of k on r is 
effectively marginalized over when the posterior distribution 
for a is calculated, so long as any sampling of radii utilizes 
the same radii at every value of a. The likelihood for a given 
power law (i.e. a given a) is thus: 



_C(a) = Y\ kr i~ a 



(15) 



where nsat is the number of satellites - i.e. the 27 compan- 
ions of M3 1 listed in Table [2] As discussed in 34.21 there 
are essentially two ways we can determine the likelihood of 
a given a. The most straightforward is to use single values 
of r,- as determined directly from the mode in the posterior 
distribution for each satellite using Eq. [7] The second, and 
arguably more robust method is to use the entire radius prob- 
ability distribution (RPD) for each satellite. In the case of this 
second approach, the likelihood for the power law determined 
for each satellite becomes a convolution of the power law with 
the satellite's RPD, so that the likelihoods of the individual 
samples are summed. The final likelihoods determined for 
each satellite can then be simply multiplied as before, giving 
a total likelihood as follows: 



£(a) = (^fr+Kr + --- + K.ii)>< 



X . . . X 



\ K '\,nsat + Kr 2,nsat + 



, . + kr 



2-a 



(16) 



nz 



kr 



2-a 



where rp is the f h sampled radius of the i th satellite, and nsam 
is the total number of samples. 

The resulting distribution achieved by implementing the 
first approach is presented in Fig. [T3](a), from which a value 
for a of 1.92*9'^ is obtained. It is interesting to note that 




Fig. 13. — Probability distributions for the slope a of a single power law used 
to model the M3 1 halo satellite distribution, given the entire set of 27 M3 1 
companions presented in Table|2] Figure (a) gives the distribution assuming 
a single best fit radius for each of the satellites as determined from the mode 
in the satellite's distance posterior distribution (as given in Column 4 of Table 
|2)- Figure (b) shows the same distribution when the entire radius probability 
distribution for each satellite is sampled 500000 times. 



this value is consistent with an isothermal satellite distribu- 
tion with uniform velocity dispersion. Replacing the indi- 
vidual best-fit radii with 500000 samples from the respective 
RPD for each satellite as per the second approach, the result is 
substantially different, as demonstrated by Fig. Q~3](b). Here 
a value for a of l^^S'^ provides the best fit to the data. 
This discrepancy is presumably a consequence of the non- 
Gaussian RPD profiles for the more poorly populated satel- 
lites, as noted in 34.21 In fact, if the 15 most Gaussian-like 
distributions are taken alone, namely Andromedas I, II, III, 
V, X, XVI, XVII, XVIII, XX, XXI, XXIII, XXIV, NGC147, 
NGC185 and M33, the results are in much closer agreement, 
with a = 1.87^y 42 with sampling and a — 2.02^^ without. 

Given the obvious asymmetry in the satellite distribution in 
Fig. [TUJ it is interesting to consider the effects of isolating 
various other satellites from the calculations. The stark asym- 
metry between the number of satellites on the near side as op- 
posed to the far side of the M3 1 tangent plane for instance (as 
had been initially reported by iMcConnachie & Irwinl [20061) 
is echoed in the respective density profiles, with an a of 



2. 37^]\^ (nosampling) recorded when only the near-side 
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satellites are considered, and that of 0.93^ ^ (nosampling) 
when instead the far-side galaxies alone are included. When 
the individual satellite RPDs are sampled, the corresponding 
values are 1.87^^ and 0.78^^ respectively. Despite the 
large uncertainties, the results clearly do not support sym- 
metry of any kind about the tangent plane. It is impor- 
tant to note however, that this asymmetry may not be phys- 
ical, but rather an effect of incompleteness in the data at the 
fainter magnitudes of t he sat ellites on the far side of M31. 
iMcConnachie & Irwinl (12006) do however observe this asym- 
metry even when only the more luminous satellites are con- 
sidered. In time, it is hoped that the nature of the data incom- 
pleteness will be better understood and effort is underway to 
determine the completeness functions for dwarf galaxy detec- 
tion in the PAndAS survey (iMartin et al.ll2012Bl) . In the mean 
time, it would seem prudent to regard the contribution to the 
density profile of the far side satellites with caution, instead 
taking the density profile measured from the near-side satel- 
lites alone as the best measurement. 

On a final note with regard to near-side-far-side asymmetry, 
it is important to realize that the uncertainty in the distance to 
M3 1 has a large effect on how many satellites will lie on ei- 
ther side of the M3 1 tangent plane, and indeed on the density 
measurement as a whole. Where the individual PPDs are sam- 
pled, this is taken into account as the M3 1 PPD is sampled for 
each measurement. Nevertheless, it is interesting to consider 
the specific (non-sampled) case where M3 1 is measured at a 
closer distance, while all best-fit satellite distances remain un- 
changed. From the M31 PPD in Fig. [121 . it can be seen that 
there is a 5% chance that M31 lies at 750 kpc or closer. If 
M31 is taken to lie at 750 kpc, Andromedas XI, XIII and XIV 
move onto the far side of the M3 1 tangent plane, going some- 
way to even out the asymmetry. However, if the distances of 
all the satellites from M3 1 are re-measured for this new M3 1 
position, the same stark contrast between the density profiles 
for the near and far sides remains and in fact grows. Using 
only those satellites on the near side of the new M3 1 tangent 
plane, an a of 2.87^q|2 is determined whereas if only those 
satellites on the far side are considered, an a of 1.22*2^2 is 
obtained. Hence it would seem unlikely that the observed 
near-side-far-side asymmetry is primarily a consequence of 
an overestimated M31 distance. 

Rec ent res earch, such as that prese nt ed by iKoch & Grebel 
(2006) and iMetz. Kroupa. & Jerienl (120071) point toward 
highly significant planar alignments of various collections of 
satellites within the M31 halo, even though as a whole, no 
such distribution is prominent. Interestingly, the former inves- 
tigation finds that it is predominantly the objects morphologi- 
cally similar to the dwarf spheroidals in their sample that can 
be constrained to a relatively thin disk, which also includes 
NGC147 and M33. While our sample is considerably larger, it 
nevertheless consists nearly entirely of such objects, so it will 
be interesting to determine what degree of symmetry may be 
found within and on either side of the best-fit plane. We intend 
to investigate this in an upcoming publication, though it must 
still be noted that outliers from the planar trend have already 
been noted in this small sample, such as Andromeda II and 
NGC185. Furthermore, other members are known not to con- 
form to the norm of M3 1 satellite dynamics, with Andromeda 
XIV for instance apparently at the escape velocity for the M3 1 



system for its determined distance (Maiewski et al. 2007). In- 
deed, it would seem that whatever model is assumed, a few 
outliers are inevitable. 

5. CONCLUSIONS 

With the ready applicability of the TRGB standard candle 
to almost any of our galactic neighbors, there can be no ques- 
tion that its role will continue to be an important one. As the 
world's premiere telescopes grow in size, so too will the ra- 
dius of the 'neighborhood' of galaxies to which the TRGB 
can be applied. Hence a technique which accurately charac- 
terizes the true probability space of the TRGB distances de- 
termined is a great asset. Indeed this quality comes to play an 
increasingly important role as more and more sparsely popu- 
lated objects are found to frequent the environs of our larger 
nearby neighbors. The differences in the results achieved in 
the previous section with and without sampling of the actual 
distance distributions illustrates this fact. 

Where in Paper I, the foundations were laid for a TRGB 
method with such desirable qualities, its full value only be- 
comes apparent when one actually employs its full Bayesian 
potential. It only requires a brief glance at figures [2] and [3] to 
see how powerful a single data-specific prior can be . Si mi- 
larly, the simple distance weighting prior outlined in {s!3.3l can 
make a poorly constrained model quite workable, as illus- 
trated in the case of Andromeda XIII. Both tools will likely 
prove very useful when the method is used further afield. 

It should also be remembered that the TRGB standard can- 
dle is in many ways, just the 'first assault.' When photometric 
data of sufficient depth is obtained, the horizontal branch can 
often pin down the distance with still greater accuracy. With a 
simple adjustment to the model luminosity function, the tech- 
niques outlined in this paper and its predecessor commute 
quite readily to implementation on the horizontal branch. 

Lastly, It must also be said that the distances presented 
herein provide an excellent opportunity to provide a new, up- 
dated analysis of the asymmetry and density of the M3 1 halo 
satellite distribution, one only touched on here. With such 
comprehensive and consistent coverage, there is great poten- 
tial in these distances to further constrain the possible evolu- 
tion and dynamical history of the M31 halo system. 
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